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Abstract 

We report two-dimensional simulations of strongly vibrated granular materials without gravity. The coefficient of 
restitution depends on the impact velocity between particles by taking into account both the viscoelastic and plastic 
deformations of particles, occurring at low and high velocities respectively. Use of this model of restitution coefficient 
leads to new unexpected behaviors. When the number of particles N is large, a loose cluster appears near the 
fixed wall, opposite the vibrating wall. The pressure exerted on the walls becomes independent of N, and linear in 
the vibration velocity V , quite as the granular temperature. The collision frequency at the vibrating wall becomes 
independent of both N and V, whereas at the fixed wall, it is linear in both N and V. These behaviors arise because 
the velocity-dependent restitution coefficient introduces a new time scale related to the collision velocity near the 
cross over from viscoelastic to plastic deformation. 
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1. Introduction 

A granular material is called a "granular gas" 
when the individual grains do not stay in contact 
with one another but rather always move separately 
through space, interacting only through dissipative 
collisions. The absence of enduring contacts between 
the particles allows granular gases to be treated by 
special numerical and theoretical methods, such as 
event driven simulations and kinetic theory. These 
methods have given rise to a large body of knowl- 
edge about dissipative granular gases LU2J ■ 

The experimental realization of granular gases is 
however more difficult, because the grains must be 
continuously supplied with energy. This is usually 
done with a vibrating plate 3J. These experiments, 
however, often result in situations quite unlike those 
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considered theoretically. Some experiments have 
then been carried out in microgravity [415] to limit 
the number of parameters in the problem in order 
to make easier the comparison with kinetic theory 
of dissipative granular gases. However, the interac- 
tion between experiments and theory has remained 
sporadic. 

In our past work [6] , we have been able to compare 
simulations with experiments by studying a gran- 
ular gas generated by placing grains in a box and 
then vibrating one of the walls to supply energy. We 
have found that a velocity-dependent restitution co- 
efficient is necessary to bring simulation and exper- 
iment into agreement. In this paper, we investigate 
a further consequence of this model that should be 
observable in vibrated granular gases in zero gravity. 

To see why the velocity-dependent restitution co- 
efficient has a radical effect in microgravity, one must 
enumerate the parameters describing the system. 
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The parameters are the following: the number of 
particles iV, the diameter of the particles d, there 
mass m, the volume of the container, the restitution 
coefficient r, the vibration amplitude A, and vibra- 
tion frequency /. Note that only one of these quan- 
tities - the inverse of the vibration frequency, 1// 
- has the dimensions of time. All other quantities 
are either dimensionless or involve length or mass. 
Thus dimensional analysis can be used to determine 
the dependence of every quantity on /. For exam- 
ple, both the average kinetic energy of the particles 
and the pressure vary as f 2 . However, this scaling is 
not in agreement with the one observed during mi- 
crogravity experiments [4] . 

There are two ways to disrupt this role of /. First, 
one could introduce gravity, bringing in a second 
time scale. The second is to introduce a velocity de- 
pendent restitution coefficient. As we show in this 
paper, this is sufficient to radically alter the behav- 
ior of the system. Specifically, the coefficient is as- 
sumed to change its character at a specific value of 
the impact velocity vq. For velocities lower than vq, 
collisions dissipate little energy, but above vq, much 
energy is dissipated. This leads to several unusual 
features that could be observed experimentally. 

2. Parameters of the simulations 

2.1. The variable coefficient of restitution 

The most important parameter in our simulations 
is the coefficient of restitution. The restitution coef- 
ficient r is the ratio between the relative normal ve- 
locities before and after impact. In contrast to most 
previous numerical studies of vibrated granular me- 
dia |7|8|9j . we let it depend on impact velocity. In 
most simulations of strongly vibrated granular me- 
dia, the coefficient of restitution is considered to be 
constant and lower than 1. 

Dissipation during collisions of metallic parti- 
cles can occur by two different mechanisms. When 
the impact velocity v is large (v > 5 m/s [TP]), 
the colliding particles deform fully plastically and 
r cc v 
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, as reported experimentally [11112113114) 
and theoretically [10112115116] . When v < 0.1 m/s 
[10j . the deformations are elastic with mainly vis- 



coelastic dissipation, and 1 — r oc v 
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as reported 



experimentally [14117118] and theoretically [17119] , 
Such velocity-dependent restitution coefficient 
models have recently shown to be important in 
numerical [20|21|22I23|24I25|26| and experimental 
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Fig. 1. The restitution coefficient r as a function of impact 
velocity v, as given in Eq. (TTJ (solid line). The dashed lines 
show vq = 0.3 m/s and tq = 0.95. Experimental points (•) 
for steel spheres were extracted from Fig.l of Ref. 1281 




[18127] studies. Applications include: granular gases 
[B], granular fluid-like properties (convection [21] , 
surface waves [22]), collective collisional processes 
[18123124] , granular compaction [23] , and planetary 
rings [25127] , 

In this paper, we use a velocity dependent restitu- 
tion coefficient r(v) and join the two regimes of dis- 
sipation (viscoelastic and plastic) together as simply 
as possible, assuming that 



(1) 



where vq = 0.3 m/s is chosen, throughout the paper, 
to be the yielding velocity for stainless steel particles 
[10128] for which r is close to 0.95 [25]. Note that 
v ~ 1 /y/p where p is the density of the particle [10] . 
We display in Fig. Q] the velocity dependent restitu- 
tion coefficient of Eq. ([T]), with ro = 0.95 and vq = 
0.3 m/s, that agrees well with experimental results 
on steel spheres from Ref. [25] . As also already noted 
by Ref. [lOj . the impact velocity to cause yield in 
metal surfaces is indeed relatively small. For metal, 
it mainly comes from the low yield stress value (Y ~ 
10 9 N/m 2 ) with respect to the elastic Young modu- 
lus (E ~ 10 11 N/m 2 ). Most impacts between metal- 
lic bodies thus involve some plastic deformation. For 
more informations on restitution coefficient mea- 
surements, see Ref. [11|12|13|14|17|18|27|28|29] , 
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Note that other laws for the velocity dependent 
restitution coefficient have been studied in the 
context of rapid granular shear flows [30] . It was 
shown that such a coefficient changes the scaling 
relation between the imposed shear rate and the 
shear stress. Specifically, when the restitution coef- 
ficient strongly decreases at high impact velocities, 
the pressure scales with the shear rate frequency 
with an exponent less than two. This finding antic- 
ipates our results, but a detailed comparison is not 
possible because Ref. [30] does not use Eq. (fTJ). 

2.2. The other simulation parameters 
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Fig. 2. Typical aspect ratio of the container for n = 5 layers 
of particles (2 mm in diameter) leading to a height of the 
bed of particles at rest, ho = 1 cm. The container height 
is h = 2.5 cm, its length L = 20 cm, and the peak-peak 
vibrational amplitude A = 2.5 cm (see text for details). 



The numerical simulation consists of an ensem- 
ble of identical hard disks of mass m«3x 10~ 5 kg 
excited vertically by a piston in a two-dimensional 
box, in the absence of gravity (see Fig. [2]). We use 
the standard event-driven simulation method, where 
collisions are assumed instantaneous and thus only 
binary collisions occur. To avoid inelastic collapse - 
an unbounded number of collisions in finite time [31j 
- collisions are made energy-conserving whenever 
very tight clusters of three particles are detected. 
Furthermore, note that using the restitution coeffi- 
cient given in Eq. ([lj prevents inelastic collapse [24] . 
For simplicity, we neglect the rotational degree of 
freedom. Collisions with the walls are treated in the 
same way as collisions between particles, except the 
walls have infinite mass. The simulation parameters 
are chosen close to the usual ones used in the ex- 
periments (see for instance Ref. [3]). The particles 
are disks d = 2 mm in diameter with stainless steel 
collision properties through vq and ro (see Fig. [T]). 
The vibrating piston at the bottom of the box has 
amplitude A = 2.5 cm, and frequencies 8 < / < 30 
Hz. The piston is nearly sinusoidally vibrated with 
a waveform made by joining two parabolas together 
6 , leading to a maximum piston velocity given by 
V = AAf in the range 0.8 < V < 3 m/s. For the pa- 
rameters used in this paper, quantities such as the 
pressure are sensitive to V, but not to the maximum 
piston acceleration [B] , so V will be used to charac- 
terize the vibration. 

The box has width L = 20 cm and horizontal pe- 
riodic boundary conditions. The number of layers of 
particles is n — Nd/L. Note that when n < 1, it is 
also the fraction of the surface covered by particles, 
so it could also be considered as an average surface 
fraction. A layer of particles, n = 1, corresponds to 
100 particles. We checked that n is an appropriate 
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Fig. 3. Snapshots of various numbers of layers n: top n = 1, 
middle: n = 5, bottom: n = 10. All snapshots taken when 
the vibrating wall (bottom) reaches its lowest point. 



way to measure the number of particles by also run- 
ning simulations at L = 10 cm and L = 40 cm. None 
of this paper's results depend significantly on L. The 
height h of the box depends on the number of parti- 
cles in order to have a constant difference h — ho = 
1.5 cm, where ho is the height of the bed of particles 
at rest, h being defined from the piston at its highest 
position (see Fig. [5]). h — ho is keep constant during 
most of the simulations (except when notified). All 
the simulations performed here are without gravity 
(except when notified) . 
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3. Results of simulations 

3.1. Snapshots 

The simplest way to display the results of the sim- 
ulation is simply to show the positions of the parti- 
cles. This is done in Fig. [3] for three simulations at 
different particle numbers. In all three panels, the 
wall shown at the top is fixed, and the bottom, vi- 
brating, wall is is shown at its lowest position. The 
periodic boundary conditions are shown by dotted 
lines at the sides of each picture. Note that through 
out the paper we will use words such as "upper" , 
"lower" , "horizontal" , and "vertical" as suggested 
by the orientation of these figures, although there is 
no gravity. 

In the upper panel (n = 1), the system merits the 
name of "granular gas" : the particles are well sep- 
arated. For the middle (n = 5) and lower (n = 10) 
panels, the situation has changed. A dense cluster 
forms against the stationary wall. As more particles 
are added, this cluster simply becomes thicker. Since 
the distance between the vibrating and stationary 
walls grows with the number of particles, this pro- 
cess can continue indefinitely. 

3.2. The pressure 

Next we concentrate on the pressure on the up- 
per wall (opposite the vibrating piston) since this is 
the quantity most accessible to experiments. Here, 
the pressure is defined as the force per unit length 
that the particles exert on this wall, or alternatively 
the flux of momentum, per unit time and length, 
through this wall. Since the collisions are instanta- 
neously, a precise temporal resolution would yield a 
series of delta- function peaks. We average the pres- 
sure over many (6400) cycles to obtain a stable av- 
erage. 

We present in Fig. |4] the dependence of P on the 
piston velocity, V , and on the number of layers n. 
This figure can be divided into two parts. The dom- 
inant feature is the pressure peak observed near n w 
1. A similar peak appears in the presence of grav- 
ity [6j . It is related to the increase of pressure with 
n up to n < 1 since interparticle collisions are rare 
and most of the particles are in vertical ballistic mo- 
tion between the piston and the lid [5] . On the other 
hand, for n > 3, the pressure is approximately in- 
dependent of the number of particles (see Fig. |4|) , 
and proportional to the piston velocity as shown by 
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Fig. 4. Time averaged pressure P on the top of the cell as a 
function of particle layer, n, for various vibration velocities: 
V = 0.8 to 3 m/s with steps of 0.2 m/s (from lowest to 
uppermost curve). Inset: P as a function of V for n = 10. 

the inset of Fig. [3] The reason for this is discussed 
below. Under gravity, the situation is quite different 
[B] . Adding particles causes more frequent interpar- 
ticle collision, and the energy dissipation to increase 
and thus the pressure to decrease [6] . At large values 
of n, resonances also appear under gravity. 

In Fig. El we examine how changing h — ho affects 
the pressure. First, we note that the independence 
of the pressure on n is not confined to a special value 
of h— ho, but holds for all values, except the smallest 
(h — ho = 0.5 cm). The pressure decreases as the 
height increases. Examining the dependence of the 
pressure on h — ho shows that P oc (h — ho) as 
displayed in the inset in Fig. [5] 

3.3. Scaling relations for global guantities 

We would now like to present the information pre- 
sented in the previous sections in a more compact 
way, and also consider other global quantities. In ad- 
dition to the pressure, we will examine the granular 
temperature T, or mean kinetic energy per particle, 
and the collision frequency C up of particles with the 
upper wall. From these quantities, one can deduce 
the average impulse AI of a particle-wall collision, 
AI up = P/C U p. 

The relation between the global quantities and the 
vibration velocity V is reasonably well described by 
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Fig. 5. Pressure as a function of n, for various heights of the 
cell: 0.5< h — ho < 5 cm with steps of 0.5 cm (from lowest 
to uppermost curve). Inset: P as a function of (h — hg)~ 1 
for n = 10. 



power-laws 




y0(n) 
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where the exponents 9(n) depend on the number of 
layers [B]. These exponents are shown in Fig.[B^, for 
a constant restitution coefficient, and in Fig. [B}d for 
a velocity-dependent restitution coefficient. 

Each point in these figures was obtained by fixing 
n and performing eleven simulations, varying V from 
1 to 3 m/s. The amplitude A is kept constant, and 
thus the vibration frequency is proportional to V, 
as explained in Sec. 12.21 Then log A (where X is the 
quantity being considered) is plotted against log V. 
The resulting curve is always nearly a straight line. 
(Specifically |log A obsorvcd - logA fitt e d | < 0.1 for all 
points). The slope of this line gives the power-law 
exponent 9. 

For the case of constant coefficient of restitution, 
?' = 0.95 (see Fig.|6^), the scaling exponents are in- 
dependent of the number of particles: C up ~ V 1 , P 
and T ~ V 2 . As said above, this is precisely what is 
to be expected from dimensional analysis, since the 
vibration fixes the only time scale in the problem. 
On the other hand, when r = r(v), a more compli- 
cated behavior is observed (see Fig. [5b). When n 
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Fig. 6. Exponents as a function of n for simulations with 
(a) r = 0.95 or (b) r = r(v). The granular temperature, T 
(0), collision frequency, C up (*), and pressure, P (o), scale 
like V e(n '>. 

is small, the exponents approach those of the pre- 
vious case. However, when the number of particles 
becomes large, all exponents approach unity. The 
reason for this is discussed below. 



4. Anomalous scalings in the dense regime 



In Sects . [3721 and 13731 we have shown that, with a 
velocity dependent restitution coefficient, the pres- 
sure of a dense granular gas without gravity obeys 
the simple non-extensive relation 



h — ho 

Let us now try to explain below this scaling. 



(3) 
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4.1. Collision frequencies at the walls 



In Figs.[3J>c, we observe that the majority of par- 
ticles remain in a loose cluster pushed against the 
stationary wall, opposite the piston. Only those par- 
ticles that "evaporate" from the cluster are struck 
by the piston. The flux of evaporating particles can 
be estimated from the rate Ci ow of collisions between 
the piston and the particles. This collision rate has a 
very curious behavior, as observed in Fig. [7] C ow is 
roughly independent of n when n > 3 (see Fig.[7K). 
This behavior holds for other values of h — ho (see 
Fig. [TJd). The inset of Fig. [7t> shows that Ci ow oc 
(h — ho) -1 . Moreover, the inset of Fig. [7^ shows 
that Ciow does not significantly depend on V (see 
the scale on the y-axes), and can be approximately 
considered as being independent of V. Thus, at high 
enough density, the collision frequency on the vi- 
brating wall is found to be 



Clow oc 



N°V° 
h — ho 



(4) 



The rate C up of collisions between the particles 
and the fixed wall is displayed in Fig. [5] It has a 
quite different behavior from Ci ow • In the dilute limit 
(n < 2), C up increases more slowly than n as al- 
ready observed in microgravity experiments [5] and 
simulations [32j . In the dense regime, when n > 3, 
C up is proportional to both the number of layers n 
(see Fig. [8]), and the piston velocity V (see inset of 
Fig- IH]) ■ Thus, at high enough density, the collision 
frequency on the fixed wall is found to be 



C up oc . 



(5) 



4.2. Explanation of the pressure scaling 

The time averaged pressure on a wall is the time 
averaged momentum flux divided by the area of a 
wall, that is 



Plow oc Cow (flow) and P up cc C up {v up } 



(6) 



and (u U p) are, respectively, the mean 
particle velocities close to the piston and close to 
the stationary wall. Since momentum is conserved, 
the flux of momentum entering the system at the 
piston, Cow (flow), must have the same value that 
the one leaving the system through the stationary 



wall, C UP 



, Therefore, the pressure on both sides 



(a) h — ho = 1.5 cm 
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(b) V = 2 m/s 
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Fig. 7. The particle-piston collision rate Ci ow as a function 
of particle layer, n. (a) At fixed h — ho, for various vibration 
velocities V = 1 to 3 m/s with steps of 0.2 m/s (from lowest 
to uppermost curve). Inset: Ci ow vs. V for n = 10. (b) At 
fixed V, for various heights h — ho = 1 to 5 cm with steps 
of 0.5 cm (from lowest to uppermost curve). Inset: Ci ow vs. 
{h - ho)' 1 for n = 10. 



is equal and is denoted P. 
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Fig. 8. Cup the collision frequency of particles with the upper 
wall as a function of n for for vibration velocities V = 0.8 m/s 
(lower curve) to 3 m/s (upper curve) with steps of 0.2 m/s. 
Inset: C up vs V for n = 10. 

Figure[3t shows that few particles are evaporated 
from the cluster, and are close to the piston. The 
evaporated particles from the cluster have the typi- 
cal velocity of the particles within the cluster, that 
is vq (see below), and thus does not depend on the 
piston velocity. Moreover, there is no reason that 
the number of evaporated particles depends on N. 
Therefore, one can expected that C\ ow oc N°V° 
which is in agreement with our numerical results of 
Eq. 0. 

The scaling Ci ow oc (h — ho)' 1 probably occurs 
because a particle that evaporates from the cluster 
must travel a certain distance before it encounters 
the piston. This distance increases with h — ho and 
thus the particle's travel time also increases. During 
its voyage, the evaporated particle could be struck 
by another particle that has just encountered the 
piston. If this happens, both particles are scattered 
back into the cluster. Thus the evaporated parti- 
cle never reaches the piston. If we assume that the 
probability of an evaporated particle being scattered 
back into the cluster is independent of time, the 
number of particles reaching the piston is inversely 
proportional to h — ho- 

When these evaporated grains collide with the 
piston, they acquire an upwards velocity propor- 
tional to V. Thus, the mean particle velocity close 
to the piston is (v\ ow ) oc V. Therefore, using these 



two above results, one have P\ ow oc Ci ow (i>i ow ) oc 
N^V 1 / {h — ho) in agreement with our numerical re- 
sults of Eq. ®. 

Similarly, the particles close to the fixed wall are 
within a cluster (see Fig.[3h) . Due to their numerous 
dissipative collisions within the cluster, these parti- 
cles move little, even less that there are many par- 
ticles within the cluster. Their mean velocity (v up ) 
is thus fixed by the dissipation within the cluster 
(thus by r(v) via vq), and by the number of parti- 
cles within the cluster (that is by N). Thus, close 
to the upper wall, one expect (v np ) oc vo/N. When 
the cluster is pushed against the upper wall, each 
layer contributes a fixed number of collisions. Thus 
the number of collisions per cycle is proportional to 
N. The collision rate is also proportional to V, be- 
cause the number of cycles per unit time grows lin- 
early with V. Thus, one have C up oc TV 1 !/ 1 which is 
in agreement with our numerical results of Eq. ([5]) . 
Therefore, using these two above results, one expect 
P up oc C U p{v up ) oc TV ]/ 1 in agreement with our nu- 
merical results of Eq. Q . 

5. Is the granular temperature relevant for 
dense granular gases? 

In this section, we verified that the notion of gran- 
ular temperature (mean kinetic energy per parti- 
cle) is relevant in our dense system. Generally, for 
an homogeneous dilute granular gas, the pressure is 
proportional to the temperature. We wonder if the 
pressure scaling with the piston velocity, P oc V, 
can be extended to the granular temperature for our 
dense system: T oc VI Although our system is not 
spatially homogeneous and not stationary during a 
vibration cycle, we will see that the linear depen- 
dence of the granular temperature on V is however 
meaningful as described below. 

5.1. Temporal distribution of energy 

We now examine the behavior of the simulations 
more closely, focusing on the variation of kinetic 
energy within one vibration cycle. We define the 
"phase" of the vibration to be a number between 
and 1 that gives the position of the vibrating wall. 
When the wall is at its lowest position, <f> = or <f> — 
1. When it has reached its highest position, (j> = 1/2. 
When it is halfway between its highest and lowest 
position, <j> = 1/4 if it is ascending, <j> — 3/4 if it is 
descending. 
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Fig. 9. Total kinetic energy of the particles as a function of 
vibration phase <j> for three different velocities V = 1, 2 and 
3 m/s (from lower to upper) at constant number of particles 
n = 5. The wall is at its lowest at cj> = and <j> = 1, and 
reaches its highest point at <j> = 0.5. Data from 100 cycles 
were averaged to obtain these curves. Inset: Total kinetic 
energy rescaled by V as a function of rf>. 

Fig. [9] shows the total kinetic energy of the par- 
ticles as a function of the phase <j>, for three differ- 
ent piston velocities, all with the velocity dependent 
restitution coefficient. Note that the kinetic energy 
varies varies by a factor of about six for each V. 
The maximum occurs around 4> « 0.3, just after 
the vibrating wall has attained its maximum veloc- 
ity. Considering the strong variations of kinetic en- 
ergy during one cycle, one might question whether 
the granular temperature T of the system were well- 
defined, or whether the law T oc P cx V (see Eq. [3]) 
is meaningful. 

The law is meaningful, because the curves of Fig. [5] 
lie on one another if rescaled with V as shown in 
the inset of Fig. [9l Thus, the granular temperature 
as a function of phase has the form T(cj>) = g(4>)V 1 , 
where g is a function describing the shape of the 
curves in Fig. [9l If one measures T at constant </> 
the same scaling law will be observed, independent 
of <j>. Note that to obtain the scaling exponents in 
Fig. El the granular temperature was measured 20 
times per cycle, and then averaged. 

5.2. Spatial distribution of energy 

Under gravity, the altitude dependence of the den- 
sity is usually measured (in order to extract the 
granular temperature). Far enough from the piston, 
the density decreases exponentially with altitude. 
A dense upper region supported on a fluidized low- 



Fig. 10. Kinetic energy as a function of altitude y at 3 
different phases <j> of the vibrating cycle. Each panel contains 
three curves corresponding to the three different V in Fig. [9] 
Data shown are for V = lm/s (lower curve), V = 2 m/s 
(middle curve) , and V = 3 m/s (upper curve) , at constant 
number of layers n = 5. 

density region near the vibrating piston is also re- 
ported experimentally [3j , numerically |33] and pre- 
dicted theoretically [34J. 

Without gravity, the spatial-dependence of the en- 
ergy is examined in Fig. 1101 These graphs were ob- 
tained by dividing the simulation domain into strips 
of height 2 mm, and then calculating the kinetic en- 
ergy present in each strip. Since the particles also 
have a diameter of 2 mm, each one will overlap two 
different strips. A fraction of the particle's kinetic 
energy is assigned to these two strips, in proportion 
to the area of the particle located in each strip. This 
procedure is carried out for all particles at fixed val- 
ues of the phase <f>, and the results averaged over 100 
cycles. 

The top panel shows the energy when the wall is 
at its lowest point (<fi = 0). At this phase, the energy 
is very low (note that the scales on the y-axes of the 
three panels are all different), and its distribution 
resembles that of the density. In the second panel, 
the wall is now just past its maximum velocity. There 
is a peak near y w 22 mm, due to the kinetic energy 
just injected by the wall. This kinetic energy is ten 
times larger than the kinetic energy found in the 
cluster, in spite of the small density. In the last panel, 
the wall has begun to move downward. Note that in 
all panels, the kinetic energy at any point is roughly 
proportional to V. Thus if one measures the kinetic 
energy density at fixed <f> and y while varying V, one 
will observe a linear dependence on V. 

Figs. 151 and ITTJ1 present a fairly complete descrip- 
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Fig. 11. Division of kinetic energy into mean motion ( "hy- 
drodynamic" ) and disorganized parts ( "thermal" ) as a func- 
tion of altitude at 3 different phases cf> of the vibrating cycle. 
(• • ■) Thermal kinetic energy (data for V = 3 from Fig. HOC . 
(— ) Hydrodynamic kinetic energy, obtained by averaging the 
velocities of all particles found in each strip. 

tion of how energy flows through the system. Energy 
is injected for 0.2 < 4> < 0.4, as the vibrating wall 
moves upward, colliding with some particles that 
have escaped from the cluster found near the upper 
wall. These particles collide with this cluster at (j> w 
0.5, exciting this cluster. Then the energy decays, 
so that by the time the wall begins moving upward 
again, most of the energy has been dissipated. 



5.3. Hydrodynamic and thermal kinetic energy 

In this section we consider the ratio of hydrody- 
namic to total kinetic energy. Here, we use the word 
"hydrodynamic" as in the context of granular ki- 
netic theory. It does not refer to any fluid moving 
among the grains, but concerns the decomposition 
of each grain's velocity into an average "hydrody- 
namic" and a remaining "thermal" component. The 
hydrodynamic velocity of a grain is found by aver- 
aging the velocities of all nearby grains. See Refs. 
35 36J for a discussion on the distinction between 
these two energies in a granular medium. 

The fraction of energy contained in the hydrody- 
namic velocities measures the organization of the 
flow. If it is close to one, then all the grains have 
nearly the same velocity. If it is small, the granular 
medium is in a gas-like state, where the randomized 
motions of the particles dominate. The use of the 
terms "granular gas" and "granular temperature" 
imply that the granular medium is in a state like that 
of a usual gas: that the "thermal" velocities of the 
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grains are much larger than the hydrodynamic ones. 
But is this really the case? One could easily imag- 
ine a situation where a cluster of particles bounces 
back and forth between the two walls, without much 
relative motion between neighboring grains. 

To see what situation applies to our simulations, 
we return to the data used to produce Fig. [TQl The 
average velocity of the particles in each strip can 
be calculated, and the kinetic energy related to this 
motion can be compared to the total. The results 
are shown in Fig.[TT] The results show that near the 
piston, most of the energy is in the mean motion 
of the particles, even when the piston is descending 
((f) — 2/3). Thus the motion of the piston is super- 
sonic. Near the fixed wall, where most of the parti- 
cles are located, however, most of the energy is ther- 
mal. One can therefore conclude that in the clus- 
ter near the wall, gas-like conditions do prevail, i.e., 
most of the motion is thermal. Obviously, our dense 
granular gas differs from an usual gas in many other 
respects, such as the formation of cluster near the 
wall, and to the anomalous scalings reported here. 



6. Conclusions 

We report simulations of two-dimensional dense 
granular gas without gravity vibrated by a piston. 
The restitution coefficient used here depends on the 
relative velocity of particles. This allows to simulate 
a dissipative granular gas in a much more realis- 
tic way than using a constant restitution coefficient. 
This model of velocity dependent restitution coeffi- 
cient is indeed in good agreement with experiments 
[Illl2ll3ll4ll7ll8l27l28j . At high enough density, 
we observe a loose cluster near the wall opposite to 
the vibrating one. This leads to unexpected scalings: 
the pressure, P, and the granular temperature, T, 
scale linearly with the piston velocity V. The colli- 
sion frequency at the fixed wall and at the vibrating 
one scales respectively, as iV 1 !/ 1 and N°V°, where N 
is the number of particles. We emphasize that these 
scalings can only be reproduced with this velocity 
dependent restitution coefficient. If one uses a con- 
stant restitution coefficient (as in most of previous 
simulations of granular gases) , one obtains FocToc 
V 2 without gravity, no matter the constant value 
of the restitution coefficient. However, this V 2 scal- 
ing is not in agreement with the one reported dur- 
ing microgravity experiments in a dilute regime [1] . 
Simulations of a dilute granular gas with velocity 
dependent restitution coefficient yield a scaling in 
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agreement with this experimental one [5] . 

One difference between our simulations and the 
microgravity experiments on granular gases is that 
it is common to shake the whole container filled with 
particles in the experiments [4] . One experiment has 
been recently performed by agitating dilute parti- 
cles with a piston in low gravity [5]. The anoma- 
lous scalings, reported here numerically in a dense 
regime, may thus be observable in such micrograv- 
ity experiments with many more particles. 
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